A supercircle description of universal three-body states in two dimensions 
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Bound states of asymmetric three-body systems confined to two dimensions are currently un- 
known. In the universal regime, two energy ratios and two mass ratios provide complete knowledge 
of the three-body energy measured in units of one two-body energy. We compute the three-body 
energy for general systems using numerical momentum-space techniques. The lowest number of 
stable bound states is produced when one mass is larger than two similar masses. We focus on 
selected asymmetric systems of interest in cold atom physics. The scaled three-body energy and 
the two scaled two-body energies are related through an equation for a supercircle whose radius in- 
creases almost linearly with three-body energy. The exponents exhibit an increasing behavior with 
three-body energy. The mass dependence is highly non-trivial. Based on our numerical findings, we 
give a simple relation that predicts the universal three-body energy. 

PACS numbers: 03.65.Ge, 21.45.-v, 36.40.-c, 67.85.-d 



I. INTRODUCTION. 

The spatial dimensions are crucial for the behavior of 
quantum systems. The reason for this can be under- 
stood by considering the kinetic energy only, and observ- 
ing that it behaves vastly different in two dimensions as 
compared to one and three. In fact, the presence of a 
negative centrifugal barrier for zero angular momentum 
in two dimensions implies that an infinitesimal attractive 
potential produces a universal bound state [l|,@|- This in- 
teresting non-trivial role of the spatial dimensionality can 
be isolated and studied in detail in the so-called universal 
limit where the effect of interactions can be described by 
one or two model-independent quantities and the detailed 
structure of two-body potentials at short-range thus be- 
comes irrelevant. 

The implementation of Feshbach resonances in cold 
atomic gases [|| allows experimenters to tune their sys- 
tems into the universal regime and study properties that 
are independent of the particularities of the atoms and 
molecules. One of the great successes of this program is 
the observation of three-body Efimov states [J, Q which 
had been proposed but never found in nuclear physics 
@. The presence of resonances also allows the study of 
the strongly correlated unitary limit, where the physics 
is universal and governed by a single so-called contact 
parameter Q that is measurable [1, [t| . Currently there 
is an experimental push to reduce the ultracold systems 
from three to two dimensions and study the universal be- 
havior [Iol - fl2j . The highly desired theoretical knowledge 
of strongly correlated gases is more sparse than in the 
three-dimensional case. Recently, the contact has been 
derived using different approaches [l3rll6| . The question 
of few-body states in equal mass two-dimensional sys- 
tems has been addressed in a number of works. Bound 
states of three- [171422T ] , four- [2l|, [23| , and larger particle 



number droplets [2514271 
bination observables 124 



and scattering [2l[ and recom- 
have been addressed. 



The general three-body bound state spectrum in the 
universal limit in two dimensions has not previously been 
investigated. The limiting case of three identical bosons 
is known to have exactly one two-body and two three- 
body bound states in the universal limit fl7H2rj( . and 
all ratios of energies and radii are universal constants. 
This is in marked contrast to three dimensions where 
the Efimov effect produces an infinite ladder of states 
when the scattering length diverges [28j. While the re- 
sults in three dimensions can be obtained through ana- 
lytical means, the spectrum in the two-dimensional case 
has been obtained through numerical methods only, and 
no analytical insights have been found so far. Further- 
more, the problem of three particles with different masses 
and short-range interactions in two dimensions has not 
been addressed before. The combination of three dif- 
ferent species of atoms was recently reported Al- 
though that experiment is presently in three dimensions, 
it should be easily extended to a two-dimensional setup. 

The universal limit can be considered in a natural way 
through zero-range interactions that leave only a single 
parameter in the potential to be related to the scattering 
length or two-body bound state energy. Here we nu- 
merically solve the momentum-space integral equations 
in two dimensions [30]. In three dimensions similar for- 
malism has been discussed within the language of effec- 
tive field theory [3l|. We extract universal features of 
completely asymmetric three-body systems and provide 
simple parametrizations of energies as functions of the 
two-body properties in terms of Lame curves 32] . This 
provides a uniform paramctrization of all aspects of the 
three-body bound state spectrum in two dimensions that 
can aid the search for a (semi-)analytical understanding 
of these intruiging systems. 
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relations: 

F n (e bc ,e ac ,m a ,mp) = F n (e ac , e bc , mp, m a ) 

f jp (A. hi — IM^ - , T? 
\e bc e bc m a m a ) 



(2) 
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The energy and mass scaling leaves us with four non- 
trivial parameters. Using the symmetry relations in 
Eq. @ we can restrict the investigations of F n to smaller 
regions of this four-parameter space as indicated by the 
dashed symmetry lines in Fig. [T] 



FIG. 1: Phase diagram of the maximum number of three-body 
bound states as functions of two mass ratios, m a and mp. The 
three two-body energies are equal, E a b — E bc — E ac . The 
roman numerals indicated the maximum number of bound 
states in each region. The circular ( 6 Li 40 K 87 Rb) and diamond 
( 40 K 87 Rb 133 Cs) points are the realistic cases studied in Figs. [3] 
and |4] The two sets of points are related by the symmetries 
in Eq. ©. 



II. PROCEDURE AND DEFINITIONS. 



We consider three particles, a, b, and c of mass ra a , m b , 
and m c with pairwise contact interactions which produce 
two-body bound states of energies E ab , E ac , and E bc . 
The wave equation in momentum space is a set of three 
coupled integro-differential equations [13, H(J HU ■ Here 
we provide a full survey of the three-body energies, E3, 
in the most general setup. Initially £3 is a function of 
six parameters. However, the use of two-body energies 
instead of interaction strengths implies that only mass 
and energy ratios enter the equations. This means that 
63 = E 3 /E ab divided by one of the two-body energies can 
be expressed as a function of four dimensionless param- 
eters, i.e. 



III. SURVEY OF MASS DEPENDENCE. 

The mass dependence of the number of three-body 
bound states for a system with E ab = E ac = E bc is 
shown. In the central region around equal masses we 
have the two bound states [T3|. This region, labeled II, 
extends in three directions corresponding to one heavy 
and two rather similar light particles, that is either 
m a ~ 1, mp ~ 1, or m a ~ rap < 1. Moving away 
from these regions in Fig. [TJ the number of stable bound 
states increases in all directions. As an example, con- 
sider m a = 10 and vary mp from small to large value in 
Fig. [TJ Perhaps surprisingly, along this line the number 
of bound states decreases to a minimum of two and sub- 
sequently increases again. The reason is that the disap- 
pearing states merge into the two-body continuum. The 
same behavior is found in three dimensions for the dis- 
appearance of the infinitely many Efimov states when 
the attractive strength is increased 34j . Variation of the 
two-body energies from all being equal leads to distortion 
of the figure but the main structure remains. The cen- 
tral region still has the smallest number of stable bound 
states. Any variation in the two-body energies make re- 
gion II bigger, pushing the other lines away from the cen- 
ter. Thus, Fig. [JJ shows the maximum number of stable 
bound states for one system described by (m a , mp). 



A. Three-body energies for given masses. 



/ E bc E ac m b m c \ _ 
\E ab E ab m a m a J 

(1) 



where 63 is the scaled three-body energy, m a = m b /m a , 
mp = m c /m ai e bc = E bc /E ab , and e ac = E ac /E ab . The 
universal functions, F n , are labeled by the subscript n to 
distinguish between ground, n = 0, and excited states, 
n > 0. By interchange of particle labels we find that 
all the universal functions, F n , must obey the symmetry 



Realistic scenarios correspond to given particles (atoms 
or molecules) with known masses. In contrast, the inter- 
actions are variable through Feshbach resonances . We 
therefore assume masses corresponding to alkali atoms 
87 Rb, 40 K, and 6 Li. The two ratios of two-body ener- 
gies are left as variables where each set uniquely speci- 
fies the three-body energies of ground and possibly ex- 
cited states. In Fig. [21 we show a contour diagram of 
the scaled three-body energies for the two lowest stable 
bound states. The log- log plot is very convenient for 
visibility but can be deceiving. On a linear scale the 
curves of equal scaled three-body energy would look dif- 
ferent. However, the parametrizations considered below 
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FIG. 2: Contour diagrams with lines of fixed £3 values as 
function of the two-body energies e ac and e(, c . The solid and 
dashed curves are for ground and excited states, respectively. 
Here a is 87 Rb, b is 40 K, and c is 6 Li. The dotted curves show 
where the number of stable bound states change from one 
(regions most asymmetric for small energies), to two (com- 
parable size of the two-body energies), and to three (equal 
and large energies). The roman labels I, II, and III are as in 

Fig.m 



nonetheless employ the linear scale. The chosen set of 
masses only allow one, two, or three stable bound states, 
depending on the two-body energies. The corresponding 
regions are shown by dotted curves in Fig. [2] The true 
extent of the regions cannot be seen. Both region II and 
III are closed, e.g region II continues along region III up 
to energy ratios of about 10 ±5 , and the narrow region 
III is entirely embedded in region II. Other sets of mass 
ratios could open region III and allow regions inside with 
more than three stable bound states. 

In the log-log plot of Fig. [2] we can see the two-body 
energies varying by five orders of magnitude, whereas 
the scaled three-body energies for a stable system must 
be larger than all two-body energies, and in particular 
larger than unity since it is measured in units of E a b- 
The three-body energy contours connect minimum and 
maximum two-body energies, that is zero and maximum 
two-body energies for the ground state and thresholds 
boundaries for existence of the excited states. The con- 
tours appear in regular intervals with larger values for 
increasing two-body energies. Large three-body energies 
reflect small spatial extension and therefore are less in- 
teresting as it presumably is unreachable in the univer- 
sal limit. The contours pass continuously through the 
boundaries of the different regions since the ground state 
exists without knowledge of the excited states. The con- 
tours in Fig. [2] for the excited state can only appear in 
the regions with two or more stable bound states. These 
contours therefore must connect points of the boundaries 
between regions I and II. They may cross continuously 
through region III precisely as the ground state would 
cross through boundaries between regions I and II. Simi- 



FIG. 3: Ro and Ri as function of £3 for three sets of mass 
ratios corresponding to 6 Li 40 K 87 Rb and 40 K 87 Rb 133 Cs. The 
curves for the excited state show 10-Ri (10e3) in order to fit 
within the same range as the results for the ground state. 



lar contours exist within region III but we do not exhibit 
them in this narrow strip where they are allowed. The 
scaled three-body energies are often substantially larger 
than the initial two-body energies, although both arise 
from the same two-body interactions. 



IV. PARAMETRIZATION. 

The universal functions, F n , defined in Eq. (JT]) are not 
easily found analytically. However, the contour diagrams 
in Fig. [2] suggest a simple implicit dependence in terms of 
an extended Lame curve or superellipse (3^ |. Note that 
despite the log-log scale of Fig. [2J the parametrization in 
terms of Lame curves is done with the energies on a linear 
scale. The three-body energies can be written indirectly 
by supercircles, i.e. 



4" + 4; 



(3) 



where the radius, R n , and the power, t n , are functions 
of £3 and both depend on the two mass ratios. The term 
supercircle has been adopted since there is a single axis 
parameter, R n . The smallest value of £3 is unity cor- 
responding to the two-body threshold of the ab system 
used as the energy unit. 

We choose sets of masses of current interest, and 
present the fitted radius functions in Fig. [3] for both 
ground and first excited states. In 6 Li 40 K 87 Rb a 
second excited state exists whereas equal mass and 
40 K 87 Rb 133 Cs only have stable ground and first excited 
states. The parametrized results cannot be distinguished 
from the computed curves in Fig. [21 The radius functions 
turn out to be surprisingly simple, that is essentially lin- 
ear functions of 63, and essentially independent of the 
masses. For the ground state we find a slight increase 
of slope with increasing three-body energy. Average es- 
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FIG. 4: Exponents, t\ and to (inset) in the supercircle fit for 
6 Li 40 K 87 Rb and 40 K 87 Rb 133 Cs as in Fig. 

timates are 

R {e 3 ) pa 0.74e 3 - 2.5 , ^(eg) ra e 3 . (4) 

The increasing functions reflect how the contours in Fig. [2] 
move to larger two-body energies with increasing £3. 

A. Symmetric systems. 

It is very satisfactory to observe this simple linear de- 
pendence which implies that the three-body energy in- 
creases linearly with a kind of average of the two two- 
body energy ratios. It is satisfying also because the sym- 
metric system, where all particles are identical, has this 
property where two- and three-body energies are propor- 
tional in the universal limit. To approach this limit we 
assume that a = b, and Eqs. ((3]) and (j4} imply for the 
ground states that 0.74e 3 2.5 + e ac 2 1 /* c) When a = b = 
c, we get the known ratios e 3 (2.5+2 1 /*°)/0.74 = 16.52. 
For the excited state we get correspondingly 63 w 2 1 /* 1 = 
1.267. This is achieved with t w 0.30 and h » 2.93. 

A system with particles a and c identical has been stud- 
ied in Ref. [3Cj. Here the = 1 and we can find £3 from 
Eqs. ([3]) and ((4]). We get the leading order terms for £3 to 
be e oc /0.74 and e ac for large e ac , respectively for ground 
and excited states. For small e ac we find correspondingly 
that £3 approaches the constants 3.5/0.74 and 1. These 
leading order terms can be supplemented by further ex- 
pansion to any desired order. The main dependence of 
the parametrization is consistent with the numerical re- 
sults in [H. 

B. The exponents. 

The exponents, t n , are crucial to obtain the correct 
curvature of the energy contours in Fig. [5] In Fig. U we 



show the functions obtained for the same sets of masses 
as in Fig. [31 These exponents increase monotonously 
with e 3 from 0.2 (n = 0) and 1 (n = 1) at the minimum 
value of e 3 = 1. The curves bend over and eventually 
stop when the states reach a two-body threshold and be- 
come unstable. In most cases this only happens at large 
energies where the universal properties are very unlikely. 

The absolute sizes increase by about an order of mag- 
nitude from ground to first excited state. The role of 
the exponents in Eq. © is to describe the curvature of 
the contours in Fig. [5] Thus, large t n is necessary for 
strongly bending curves. This explains the difference be- 
tween ground and first excited state, but also the overall 
increase with 63. This is especially pronounced for the 
excited states which are squeezed in between boundaries 
defined by stability towards decay to bound two-body 
subsystems. 

The contours connect boundaries with small e ac and 
large the or vice versa. The end-points are when either 
£bc or e ac is close to e 3 . The slopes are from Eq. 
found to be, de ac /deb c = — (ffce/cac)*™ - • This means 
that the contours in Fig. [5] first exhibit a small negative 
slope (small the) and then a quick bend to reach a large 
negative slope (small e ac ). All contours decrease from left 
to right. The regions become narrower for higher two- 
and three-body energies, and the contours must bend 
quickly. The implication is that t n must increase with £3 
and with the number of excited states. 

The allowed regions for a specific number of bound 
states are indicated in Fig. [T] The practical examples 
selected for Figs. |3] and 0] are found in the center, and 
in regions II and III of the third quadrant. In Fig. [2] 
the boundaries for two bound states are rather widely 
separated giving rise to the smallest exponents in Fig. 2) 
Moving towards more bound states in Fig. []] would widen 
this region and the narrow strip with three bound states 
as well. At some point one more bound state appears but 
in the process all regions widened and the exponents in 
the fits consequently decrease. The other cases in Figs. [3] 
and|4]are both in region II of Fig.[T]but the point for equal 
masses is further away from the boundary to a region III 
than the point of asymmetric masses. This is a general 
explanation of the hierarchy in the sizes of the exponents. 

The behavior of the exponents is also surprisingly sim- 
ple for each set of masses. The relatively fast increase 
at small energies in Fig. @] slows down and both to and 
t\ approach constants at large energy. For the ground 
state, this can be accurately captured by 

4-1 \ ~ £ 3° + A) ,r\ 

where oq and ot\ are the mass dependent constants ap- 
proached at large energy, see Fig. 01 The parameters, 
(pOjQio) — (0.04 — 0.06,0.3 — 0.5), exhibits a small mass 
dependence, whereas (/3o, 70) — —(0.93 — 0.95), —(0.82 — 
0.87) are slightly negative but almost mass independent. 
We emphasize that stability requires £3 > 1. The value 
of to for small £3 sa 1 is then in the range of to ~ 0.2 — 0.5 
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as required to give the limiting value of £3 = 16.52. A 
similar parametrization for the exponent corresponding 
to the excited state can be found. 

V. CONCLUSIONS AND PERSPECTIVES. 

The general spectrum for three interacting particles in 
the universal regime in two dimensions is currently un- 
known. In the simplest case of three identical bosons, 
there is a lack of analytical results and one must re- 
sort to a numerical treatment at the moment. We nu- 
merically compute the universal three-body energies for 
short-range interactions, different masses and interac- 
tion strengths. The number of bound states varies from 
one and up depending on mass ratios and two-body 
subsystem energies, with symmetric system having the 
fewest and the most for two heavy and a light particle. 
The three-body ground state energies can be successfully 
parametrized by universal functions that are so-called su- 



percedes (powers different from two) where the coordi- 
nates are the independent two-body energy ratios and the 
radius parameter is the three-body energy. The latter is 
an approximately linear function of the three-body en- 
ergy, independent of masses, while the powers of the co- 
ordinates are functions of both mass and three-body en- 
ergy. The simple threshold behavior for identical bosons 
where the two three-body bound states are directly pro- 
portional is only valid for relatively small three-body en- 
ergies. Whether this indicates a smaller universal region 
has to be investigated by either finite range models or ra- 
dius computations. Our results can be used as estimates 
of three-body energies and number of bound states, and 
as a measure of the deviation from the universal zero- 
range limit. 
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